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Abstract 

It has been predicted that the zero frequency density of states of YBCO in 
the superconducting phase can display interesting anisotropy effects when a 
magnetic field is applied parallel to the copper-oxide planes, due to the dia- 
magnetic response of the quasi-particles. In this paper we incorporate para- 
magnetism into the theory and show that it lessens the anisotropy and can 
even eliminate it altogether. At the same time paramagnetism also changes 
the scaling with the square root of the magnetic field first deduced by Volovik 
leading to an experimentally testable prediction. We also map out the an- 
alytic structure of the zero frequency density of states as a function of the 
diamagnetic and paramagnetic energies. At certain critical magnetic field val- 
ues we predict kinks as we vary the magnetic field. However these probably 
lie beyond currently accessible field strengths. 
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I. INTRODUCTION 



In 1993 Volovik |T|] predicted that the density of states of a (i-wave superconductor is dom- 
inated by delocahsed nodal quasi-particles and scales as the square root of the magnetic 
field. Since then there has been a lot of research on the role of nodal quasi-particles in 
determining the electronic density of states of the cuprate superconductors. Experimentally 
this is probed by looking at the specific heat as a function of magnetic field This is a 

bulk property and therefore complements STM and ARPES measurements which are surface 
probes. (On the other hand, it is a somewhat less direct probe of the electronic density of 
states.) Theoretically, groups have studied the effect of dirty superconductors @], the role of 
the copper oxide chains in affecting the specific heat and the anisotropic dependences on 
magnetic field direction p . In particular, the last paper predicted that the observed depen- 
dence of the density of states of an in plane magnetic field would depend on its orientation 
relative to the gap nodes. That paper only accounted for the diamagnetic response of the 
quasi-particles, as in Volovik's original paper. In this paper we generalise this to include the 
paramagnetic response of the quasi-particles. The quasi-particle paramagnetism has been 
considered in different contexts in P,p!0|. We focus on YBCO because only it has sufficient 
transport in the c direction to allow the usual London theory to apply so that vortices to 
form out of the CuO planes [Q. Other relevant papers that treat different aspects of the 
problem include UTTHT^. 

In the following section we present a general discussion of the theory and define the 



various terms. In section |ITl| we discuss the energy scales specific to YBCO. The two sub- 
sequent sections present the results for the physical case of relatively small magnetic field 
for the field perpendicular and parallel to the CuO planes, respectively. With small fields, a 
particularly simple approximation can be invoked leading to analytic results. In the parallel 
case, we predict a loss of anisotropy at sufficiently large magnetic fields. In section we 
present a more physical argument for the loss of anisotropy. We summarise our results in 
the conclusion. In the appendix, we allow for fields of arbitrary strength but only in the 
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situation where the field is perpendicular to the copper-oxide planes. 



II. GENERAL THEORY 



The approach we use is to consider a BCS-like formalism [|I4[ with a (i-wave gap. Since 
YBCO is a type-II superconductor, there are vortices above a critical field Hd and we will 
focus on the regime Hd <^ if ^ ifc2- This condition ensures that ^ <^ R <^ X where ^, 
R and A are the coherence length, vortex radius and penetration depth, respectively. In 
fact, Hci is so small (~ 100 Gauss |]T3[) and Hc2 so large (~ 100 Tesla) that one is always 
in this experimental range. Around each vortex there is a rotating superfiuid condensate. 
Quasiparticles are excited out of this condensate and are accordingly Doppler-shifted by the 
superfiuid velocity. The effect of this Doppler shift is to give a finite density of states, even 
at zero frequency, and we refer to this as the diamagnetic effect. In addition, the individual 
quasi-particles have a paramagnetic response to the magnetic field due to their intrinsic 
spin and we refer to this as the paramagnetic effect. In this paper, we are interested in the 
competition between these two effects, there being two important geometries. The simpler 
one is when the field is parallel to the c-axis {H \\ c) in which event we can study effects 
for all magnetic field values. The more physically interesting geometry, however, is that the 
field is parallel to the a — b plane {H \\ ah). In this case the diamagnetism results in an 
anisotropic dependence of the density of states on the direction of the applied magnetic field, 
with a greater value when the field points along the gap antinodal direction and a lesser 
value when the field points along the nodal direction. We show that the paramagnetism 
lessens this effect and at a critical magnetic field makes it disappear altogether. 

Due to the paramagnetic shift, the density of states is broken into two components which 
can be written semiclassically as 

iVf ^ « / dr /dk|V|5 [j^^il + Al±f^H^'- V'^ . (1) 

This is a simple generalisation of the formalism of ^ where we have included the para- 
magnetic term ±fiH. We describe each factor in turn. The (i-wave gap is given by 



Ak = Aq cos^(20) and = — Pc so that y^k + ^ pH is the quasi-particle excitation 
energy. (Here we have defined the quasi-particle energy relative to the chemical potential 
/ic.) The integral over k (with polar coordinates k and 0) is a standard trace integral used 
to determine the density of states. We will use the delta function to do the radial k integral. 
The integral over r is an average over one unit cell of the vortex lattice; r is the distance from 
the vortex core and (3 is the vortex winding angle. It is convenient to nondimensionalise the 
radial integral to p = r / R, where R is the inter- vortex distance 

a is a geometrical constant accounting for the mismatch between the circular vortices and 
the regular lattice they fill out and $0 = hc/2e is a fiux quantum. 

The factor V" = Vg ■ is the Doppler shift which depends on the condensate velocity Vg 
(which in turn depends on r) and on the Fermi momentum kp. We first consider H \\ c\ far 
from the vortex core (but still within the penetration depth scale) this is approximately 

V ^— sin(/? - 0) H\\ c. (3) 

P 

Eh = vpa^JiiH/A^Q is an energy scale associated with the diamagnetism (and Vp is the 
Fermi velocity.) This yields = 2A^/H |^ where Eh is measured in meV and H in Tesla. 
In practise, when integrating over the winding angle jS, we change variables to (3 — (p so 
the sinusoid is sin/3. There are two difference when H \\ ah. Firstly, the differently aligned 
vortices implies 

E ]^ 

sin /? sin(0 - a) H || ab. (4) 

There is an additional parameter a which describes the orientation of the field with respect 
to the axis directions in the ah plane. Specifically, when a = 0, the field points along the 
antinodal direction of the gap function while when a = tt/A, the field points along the nodal 
direction. The second difference is that the energy scale Eh is smaller. This is due to the 
fact that the vortices must form out of the copper oxide planes. This leads to a rescaling of 
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the effective Fermi velocity by yXab/K where Aj is the penetration depth in direction i. Eh 
is accordingly reduced by a factor of about 2.5 relative to the H \\ c case. 

For now, we shall keep V arbitrary and thereby include both field alignments with the 
same analysis. To begin, we need to find the values of k for which the argument of the delta 
function is zero. For the -\-fiH term this is given by 



^e^ + Al + f,H-\V\=0. (5) 

A solution may or may not exist depending on the values of 0, P and p. (In principle we 
should also allow for the solution with +|^|, but in practise there are no solutions to that 
equation.) We shall refer to the contribution arising from this term as Nq'^ ■* (the reason for 
this name is made clear below.) For the —fiH term the analysis is a little more complicated. 
Specifically, we now must allow for the possibility of two distinct roots of the form 



a + ^l-f^H±\V\ = 0. (6) 

Here we allow for both signs of \V\. One should not confuse the ± in (|^) with the ± in the 
choice of spin direction in (|l|). We shall call the respective contributions A^q and A^g ■* 
where the first sign refers to pH and the second to \V\. 

As stated, we use the 6 function to do the k integral. This leads to the following 
expressions for the three terms (where we normalise by recalling that for large fields, the 
total density of states should approach that of the normal state A^): 
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(7) 



The function P{x) equals a; if a; is real and positive, otherwise it is zero (this yields the ranges 
over which the respective solutions of Eqs.(|^) and (^) exist.) Our goal is to do these three 



integrals and thereby calculate the density of states and specific heat. It is convenient at this 
point to introduce the dimensionless quantities v = Eh/Aq and cr = fiH/Ao describing the 
energies associated with the diamagnetism and paramagnetism respectively. Mathematically 
we can treat them as independent quantities. In reality they are dependent since u oc \^ 
while 0" oc if so that a oc z/^, but this can be imposed later. 

We begin by assuming that the fields are small enough to use the "nodal approximation" , 
which amounts to replacing the gap function Ak by a local linear approximation in the 
neighbourhood of the gap nodes. This leads to a very simple integral. This approximation 
is valid in the limit of small u and a since in that event the range of contributing (p values 
in d^) is very small and a linear approximation is valid throughout the allowed integration 
range. It is convenient to first take H \\ c. This then sets up the small field H \\ ab 
calculation. The small field limit underlying the nodal approximation is physically relevant 
since it corresponds to currently accessible magnetic field strengths. However, it is also 
interesting to explore the structure for arbitrary fields; this is done in the appendix. 



III. ENERGY SCALES 

We now discuss the energy scales in the problem. It will be convenient to parameterise the 
diamagnetic and paramagnetic energies as follows. We define the coefficients a and b such 
that u = a\/H and a = bH and the functional relationship between v and a is a = bv^ ja^. 
Knowing the coefficients a and b gives us the full information about the thermodynamic 
response of the system. For example, the cross-over field Hf, where z/ = a is given as {a/bY- 
We begin by noting that for YBCO Aq is about 20 meV . The paramagnetic energy is 
simple since it is just which equals 0.058ii where H is measured in Tesla and the energy 
in meV (henceforth these units will be understood.) Together these imply b = 2.9x10^'^. 
The diamagnetic energy scale is trickier; to date it has been found directly from experimental 



measurements. This point is discussed in Refs. ||T6|] where it is argued that Eh is different 
for H II c and H \\ ab. In the first case. Eh ~ 2Ay/H so that a = 0.12 and a = 0.20z/^ . 
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The cross-over field is therefore He = 1.7xlO^T. This is enormous and we can confidently 
state that for all experimentally accessible fields, the diamagnetism completely dominates 
the paramagnetism. 

The H \\ ab case is more ambiguous. Based on the ratio of the penetration depths 
in the c direction and the a and b directions, the authors of Refs. [0] estimated that 
Eh ~ 0.7V^ leading to a = 0.035 so that the cross-over field is He = 147T. On the other 
hand, experiments which compared the scaling of the specific heat with \/H in the c and 
a — b directions led to the conclusion that Eh < OA^/H [|1^. Taking the upper limit of the 
inequality for concreteness, we would have a = 0.020 so that the cross-over field is He = 48T. 
While the cross-over happens at rather large fields, we will show that the loss of anisotropy 
is apparent at smaller fields which are certainly experimentally accessible. 



IV. H II c 

We first consider the field parallel to the c-axis. For this situation, all four gap nodes 
contribute identically. Also, we can use symmetry to restrict the range of (3 and thereby 
dispense with the absolute value signs around V. Using the nodal approximation to do the 
integral, we find 

-4=^^- dp d(3P{usm(3-ap) 
N n Jo Jo 

-4^^- dp / dpP{usmP + ap) 
i\ n Jo Jo 

-4^^- dp dp P{ap -usinP). (8) 

iV TT JO Jo 

(Recall that P limits the integration ranges so that the corresponding zero solution of the 
original delta function argument exists.) It is interesting to remark that there is no limitation 
on the integration range of Nq \ which contributes throughout the vortex unit cell for all 
values of a and u. The other two terms are restricted and in fact, their allowed ranges in the 
p — 13 space are complementary. There is a physical meaning to this. It says that in those 
regions of the vortex unit cell, it is energetically favourable for the quasi-particle to flip its 

7 



spin even though this costs it diamagnetic energy. In the hmit where the diamagnetism is 
dominant, we have only Nq^~\ In the limit where the paramagnetism is dominant, we have 
only A'^o"'*''' , that is all quasi-particles are spin down. In the regime where the diamagnetism 
and paramagnetism are competitive, the final result involves a nontrivial combination of 
both terms. These three terms can also be combined and expressed as 



■= ^ — dp dp (Iz/sin /? + ap\ + lap — usin P\) . 
N n Jo Jo 

From here the analysis is simple and summing all three terms we conclude 



(9) 
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(10) 



Clearly these two expressions agree when a = v. They also give the correct behaviour for 
either cr = or z/ = 0. For later purposes, it is be convenient to collectively express these as 



^^F(a,z.). 



The function F has a nice scaling property. Specifically, 



(11) 



F(a,z/) = z//(a/z/) 



(12) 



where 
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(13) 



Other than the prefactor in ([12|), the function F is given by a simple universal form. 

In Fig. |l] we show the function / plotted as a function of x. To aid comparison, we plot 
the X > 1 expression as a dashed line in the range x < 1. The change in functional form at 
X = 1 represents a very mild nonanalyticity which is not observable by eye. Analysis of the 
function /(x) shows that the first three derivatives are continuous at x = 1 and that the 
X < 1 branch suffers a (1 — x)^/^ discontinuity as it approaches x = 1 from below. 
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FIG. 1. The function / as a function of x. There is a change in functional form at a; = 1; we 
have continued the x > 1 curve down to a; = as a dashed hne for the sake of comparison. 



We now consider the more interesting situation where the magnetic field is parallel to the ab 
planes. We continue to use the nodal approximation and again assume that the 4> integral 
is dominated by the regions around the four gap nodes. We define 0„ = (2n — l)7r/4 such 
that A(0„) — and divide the integral into four domains centred on these values. In 
contrast to the previous case, there is a more complicated dependence since it appears in 
V itself. However, to the same order of approximation, we ignore the variation of when 
evaluating V, and simply use the value of appropriate to that domain. The development 
then follows the previous case but where we must consider each of the four integrals (as 
indexed by n) and also keep track of the a dependence. We also restrict the /3 range as 
before. The result is 



V. H II ab AND ANISOTROPY 
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% dp/ d/3 (-z/|sin(0„ - a) I sin /5 + cTp). 



(14) 



By inspection each of these integrals has the same structure as Eq. (||) but with a rescaled 
value for u. We then conclude 



(15) 



By symmetry, the n = 1 term equals the n = 3 term and the n = 2 term equals the n = 4 
term. 
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FIG. 2. The electronic density of states as a function of a for H = lOT and using the exper- 
imentally determined value a = 0.02. The upper curve is for the physical value b = 0.0029. The 
lower curve is for comparison and is without paramagnetism (i.e. b = 0.) 

As mentioned, one aspect of || a6 is that there can be an interesting anisotropy as we 
vary the angle a. The maximum value is for a = while the minimum is for a = n/A: 



No 
N 
No 
N 



^ F(a,u/V2 



a=0 



(F(a,0) + F(a,z/)) 



(16) 



a=7r/4 2 

We define AA^ to be the difference between these two values. There are two interesting 
limits that we can evaluate analytically: 



cr < z/ 
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^0 a>u. (17) 

When (7 = 0, the first of these is the result found in 0. We note that introducing the 
paramagnetism lessens the amount of anisotropy. This trend continues as a increases until, 
at the critical field where a = the anisotropy completely vanishes. This last result can be 
found with more direct physical arguments which we present in Section 0. 

A relevant field scale is given by where AA^ is a maximum. Analysis of (|l^ shows that 
this happens at Hm = O.OTifg- (Note that at this field, a = 0.26z/ so that the criterion cr ^ z/ 
in { \[7\) is roughly satisfied.) For the theoretical estimate of a = 0.035, we have Hm = 10 
Tesla while for the experimental value (and the one we trust more) of a = 0.02, we have 
Hm = 3.3 Tesla. It is this field scale which effectively controls when the overall behaviour 
becomes significantly affected by the paramagnetism. 

In Fig. ^ we show the angular dependence of the density of states for a fixed magnetic 
field as we increase the amount of paramagnetism. As we increase the paramagnetism, the 
density of states increases; this effect is particularly pronounced at a = 7r/4, as we expect. 
Consequently, the amount of anisotropy decreases, again in agreement with our general 
discussion. 
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FIG. 3. The electronic density of states as a function of v H (with H measure in Tesla) for 
o = 0.020. The upper curve in each case is for a = while the lower curve is for a = vr/4; the 
difference between the curves is the magnitude of anisotropy. In the left figure we took b = 0.0029, 
the right figure is for comparison without paramagnetism (i.e. b = 0.) 



Next we show the anisotropy as a function of magnetic field. In Fig. we plot the 
field dependence of the density of states using the experimental estimates of a = 0.02 and 
b = 0.0029. The upper curve for each pair is for the field orientation a = while the lower 
one is for a = n/A. Notice that they start with different slope but ultimately converge leading 
to a limited range of field strengths over which the anisotropy is present. For comparison we 
also show the results without anisotropy for which the two curves never converge. Also note 
that around the scale Em defined above, the two curves begin to converge. At around that 
scale there begins a noticeable nonlinear dependence of the density of states on ^/H, which 
is in principle measurable. The calculation presented here is incompatible with a density 
of states which is both isotropic and linear in \/H, thus providing a clear experimental 
signature. The data of |jT^ is marginal in that the error bars are large enough to not be in 
conflict with theory. However, a more precise measurement or one that went to larger fields 
would presumably be able to settle this question. Preliminary data from another group 
does seem to be consistent with a small anisotropy. 
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FIG. 4. The lower figure shows the absolute value of the anisotropy while the upper figure 
shows the relative value of the anisotropy {AN divided by the value at q = 0). The solid line is 
from using a = 0.020 while the dashed line is from using a = 0.035. 

In Fig. ^ we plot both the absolute anisotropy which initially grows with magnetic field 
as the overall scale of the diamagnetism increases with magnetic field. However at some field 
strength of order Hm it stops growing as the paramagnetic effects begin taking hold and 
finally vanishes when the paramagnetic effect dominates. We also plot the relative anistropy, 
i.e. the absolute anisotropy divided by the value at a = 0. This decreases abruptly with 
magnetic field (with a square root dependence for small field) so that even relatively small 
values of a can have a large affect on the anisotropy. We plot the results for two values of a 
showing the strong dependence of this affect on the diamagnetic energy scale. Unfortunately 
this is not a well established number, however we have the greatest confidence in a = 0.02 
which is shown as the solid curve in the figure. 

VI. HEURISTIC ARGUMENT FOR LOSS OF ANISOTROPY 

The anisotropy arises entirely from the diamagnetism, so it certainly makes sense that it is 
negligible in the limit that the paramagnetic effect is much stronger than the diamagnetic 
effect. However, it is not entirely intuitive that it should disappear completely when the 
effects are of comparable importance {i.e. when u = a.) In this section we present a 
heuristic explanation of the anisotropy. We begin by showing a picture of what happens 
when a small amount of paramagnetism is introduced to an otherwise purely diamagnetic 
problem. We then present the opposite limit of near perfect paramagnetism with a small 
amount of diamagnetism. We then bridge these two limits and in the process present a 
concrete explanation of why the anisotropy switches off precisely when u = a. 
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FIG. 5. Each "V" shows the dependence of the density of states on the argument lv. It is locally 
linear but with a cusp at zero argument. We show the results for the field aligned along the gap 
antinode (top two) and aligned along the gap node (bottom two). In each row, the left hand figure 
is for the = 37r/4 gap node and the right hand figure is for the (p = ir/A gap node. The solid line 
in each case shows the result with no paramagnetism (note that there is no contribution from the 
bottom right diagram in this case.) The contribution from the two nodes for a = outweighs the 
one contribution for a = 7r/4. We indicate the effect of adding paramagnetism with the dashed 
lines. Each solid line gets split into two contributions but the average of the two equals the result 
without paramagnetism, due to the linear behaviour, so there is no change to the net contribution. 
For the bottom right diagram this is not true and there is a net increase in the contribution. This 
lessens the anisotropy. 

In Fig. 1^ we show the dependence of the density of states as a function of frequency, in 
brief A^(ti;) oc |u;|. Although all of our calculations involve the physical frequency being zero, 
due to the Doppler and paramagnetic shifts, our final answers involve evaluating this density 
of states at various arguments. For example, in Eq. (|^), we can interpret the integrands 
\(Tp ± z/sin(0 — a) I as being the density of states N^u) evaluated at u = ap ± z/sin(0 — a). 
All of the interplay between the paramagnetism and diamagnetism involves arguments being 
near zero so that the cusp in the function plays an important role. For concreteness, we 
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take (3 = tt/2 and p = 1; we will discuss what happens elsewhere in the vortex unit cell in 
the subsequent discussion. In the purely diamagnetic limit, we set a = 0. For -f^ || c, all 
nodes contribute identically and there is nothing interesting to say. Instead we will focus 
on H II ah. In this case, the nodes at = 7r/4 and 57r/4 are identical as are the nodes 
at = 37r/4 and 77r/4, so we will only look at one of each of these pairs. Recall that 
the rescaled Doppler shift is \vsm{(j) — a)|. When a = 0, all nodes contribute identically 
as as indicated in the top two diagrams of the figure. When a = 7r/4, one pair of 

nodes contributes as v while the other does not. The dependence on the argument is linear, 
meaning a greater density of states for a = than for a = 7r/4 (since 2v j \/2 > u.) 

Imagine now that we add a small amount of paramagnetism, in other words set a to be 
a small but non-zero value. This induces a splitting in the argument at which we evaluate 
N{uj) to z/| sin(0 — a)\ ± a. As long as i/sin(0 — a) is larger in magnitude than a, this 
splitting has no effect. This is because N{u!) is locally linear so that the extra amount from 
one sign of a is exactly compensated for by the other sign of a. The only case that this does 
not apply to is the (p = n/A node when a = tt/A. In this event, both signs of a contribute 
equally and there is no cancellation. This implies an increase in the density of states for 
a = ir/A and no change for a = so that the anisotropy is lessened, in accord with Eq. (jl^. 
(We recall that this was all for the one point p = 1 and (3 = tt/2. When we average over 
the vortex unit cell, there are regions of small (3 where ap dominates the diamagnetic shift. 
This leads to corrections of order cr^ to the final, average density of states. This does not 
affect our general conclusion that adding paramagnetism lessens the anisotropy.) 

We now consider the opposite limit where the paramagnetism dominates. With no 
diamagnetism, we simply evaluate the density of states at ±a. This is clearly the same for 
all orientations and there is no anisotropy. We now add a small amount of diamagnetism, 
in other words set z/ to be a small but non-zero value. Analogously to before, we have shifts 
in the argument by the amount ±z/sin((/) — a). As before, the local linear dependence of 
N{uj) means that the net effect cancels and there is no change to the densities of states and 
hence no introduction of anisotropy. Again, recall that this was for f3 = tt/2 and p = 1. 
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The argument just presented actually breaks down for small p where the diamagnetism is 
greater than the paramagnetism and for which an altogether different mechanism preserves 
the lack of anisotropy. The diamagnetism is dominant in the region of the unit cell for 
which ap < z/| sin(0 — a)sin/5|; the area of this region clearly scales as v. The value 
of the density of states in this region scales as v (due to the local linear dependence of 
N{uj)). Therefore, the total contribution scales as z/^. This is small when the paramagnetism 
dominates. Furthermore it leads to no anisotropy since the resultant sum over the nodes 
is proportional to sin^(7r/4 — a) + sin^(37r/4 — a) = 1, which is clearly independent of a. 
The key points are the quadratic scaling in v and the rescaling of v by sin(0 — a). Taken 
together they lead to a completely isotropic a dependence. 

For cases where a and v are comparable, we continue to use the arguments of the above 
paragraph. As long as z/ < a, the scaling with v'^ is perfect (as is also apparent in Eq. ([T0|)) 
and there is no anisotropy. As soon as z/ > o", the scaling breaks down. This is because 
the curve ap = |z/sin((/) — a) sin/5| exceeds the maximum value of p = 1 and the area of the 
vortex unit cell in which the diamagnetism dominates begins to saturate. The contribution 
of this part of the vortex unit cell then increases more slowly than z/^ and a nontrivial a 
dependence becomes possible. Nevertheless, the scaling is still approximately maintained 
even for u > a so that the anisotropy turns on very slowly. (This is related to the very 
smooth behaviour of the function f{x) near x = 1 shown in Fig. |l|.) Thus even for a quite 
a bit smaller than u, the anisotropy is strongly suppressed. 

We remark that we made extensive use of the local linear behaviour of the density of 
states. This linear behaviour arises from using the nodal approximation for small fields. A 
more general discussion allowing for quadratic and higher powers will not affect the general 
conclusions but presumably would lead to a small amount of anisotropy for a > u. 
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VII. CONCLUSION 



In this paper we have discussed the combined effect of paramagnetism and diamagnetism on 
the zero frequency quasiparticle density of states in the presence of a magnetic field. When 
the field is parallel to the c-axis, there is no interesting anisotropic effect. Furthermore, 
for the cuprates, the energy scales are such that the paramagnetism is negligible for any 
physically realisable magnetic field strength. 

When the field is in the ah planes, the diamagnetism is less important — this is related to 
the fact that it is more difficult to form vortices out of the copper-oxide planes. In fact, we 
expect that it is the paramagnetism which is dominant for all materials but YBCO. YBCO 
is special since the c— axis transport, necessary for the formation of vortices out of the plane, 
is small but not negligeable. For this material, we expect that the paramagnetic and dia- 
magnetic effects to be of comparable importance for magnetic fields of experimental interest. 
One important issue is the anisotropy as we change the angle between the applied field di- 
rection and the gap nodes. Even a small amount of paramagnetism strongly suppresses the 
expected anisotropic response, providing a possible explanation of why this effect has not 
been observed in experiments. However, the same field scale where the paramagnetism be- 
gins to strongly suppress the anisotropy also controls where the linear dependence of on 
\/H begins to fail, as is apparent in Fig. ^ Therefore, if this is the explanation, one should 
also observe nonlinearities in such experimental curves. The experiments of fl^ seem to 
indicate both no anisotropy and linear behaviour in \/H , in contradiction with our conclu- 
sions. However, the error bars are large enough to still be consistent with our conclusions. 
More precise experiments or ones using larger fields are necessary to resolve this question. 
Preliminary data from |^ indicates that there is an observable anisotropy. 

Research supported in part by the Natural Sciences and Engineering Research Council 
(NSERC) and by the Canadian Institute for Advanced Research (CIAR). 
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APPENDIX A: GENERAL RESULTS FOR H 11 c 



In this appendix we generalise the results found using the nodal approximation by not 
considering either a or to be small but limiting ourselves to the somewhat simpler H \\ c 
case. Consider equations (|^). For a given choice of p and /5, there is always some finite 

range of such that a solution to the ( ) equation exists; this is not true for the (H — ) 

and ( — h) equations. The resulting integrals correspond to the definition of the elliptic 
function and we conclude 



, , d(3dppG{a) 

IS TC 



dpdpp G{b) 



-(+-) 4 I- 

IS TT^ Jv>aH 



where we have defined the function 



G{x) = K{x) X < 1 

^-K (-\ x>l. (A2) 



X \x 



K{x) is the complete elliptic integral of the first kind and the parameters are 

P 1. P . P 



(A3) 



z/sin/5 + crp crp — usmP i/sin/3 — crp 

There is a logarithmic singularity when the argument of G equals unity but this is integrable. 
The integration domain of the first equation in is < p < 1 and 0</3<7r/2so that 
y = |V^| > 0. The domain is the same for the other two except for the added constraint on 
V -pH. 

It is difficult to go further in an explicit analytic determination of the integrals. However, 
we can say a few things on a qualitative level. There are logarithmic singularities whenever 
one of the quantities a, 6 or c is unity. These possibilities correspond to curves of singularities 
in the p — P space. Also, when b (or, equivalently, c) equals infinity, there is a curve of zeros 
in the p — (i space. As v and a are varied, singular things can happen to these special 



curves. For example, in order to have a logarithmic singularity in A^q ■> we require that 
a = 1 which in turn requires 

p = sin (A4) 

1 — (7 

Clearly if cr > 1 there is no solution to this condition so that as a crosses through unity 
we will notice the abrupt disappearance of the curve of logarithmic singularities in the 
integrand of Nq For Nq~~^\ it is precisely the converse. There are only logarithmic 
singularities if cr > 1 and none if u < 1. Therefore, as a crosses through unity, we expect 
some nonanalyticity in the density of states. 

Another possibihty for A^q ^ is that the curve of singularities crosses through p — I 
and out of the allowed domain of the p integral. By the previous equation, this will happen 
when v — 1 — a. The analogous event happens for Nq when v — a — 1 and for Nq'^ ■* 
when u — 1 + a. The final possibility is that the curve of zero values for Nq^~^ and Nq"^^ 
crosses through p = I and out of the allowed domain of the p integral. The curve of zero 
values (when b and c are infinite) requires p = i/ sin /3/cr so this event occurs when u = a. 
We already encountered the a — u criterion in our small argument analysis and is shown in 
Fig. 1. The others all require that either u or a (or both) be of the order of unity and so 
were not captured in our small argument analysis. 
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FIG. 6. Top: A plot of the different lines in the v — a space along which we expect nonanalyt- 
icities in the density of states. The two parabolas represent possible traversals through this space 
as we increase the magnetic field. Bottom: the resultant dependence of the density of states for 
the two parabolas of the upper plot. The horizontal axis is in arbitrary units. 



A plot of the different domains in the v — a space is shown in Fig. ^. Each of the solid 
lines represents some border along which one of the three terms experiences a nonanalyticity. 
As we change the magnetic field, we follow a parabola in this parameter space and thereby 
intersect various of these lines. The precise sequence of intersections depends on the energy 
scales associated with the paramagnetic and diamagnetic terms. There are five possible 
sequences depending on the sharpness of the parabola; we show two of them. For a very 
sharp parabola, the paramagnetic energy is dominant except at very small fields and the 
result is similar to a purely paramagnetic calculation. The major differences are 1) an initial 
y/H dependence for small enough fields and 2) the logarithmic singularity at cr = 1, which 
would be present for pure paramagnetism, is softened by the diamagnetism. 

In a situation where the diamagnetism is dominant for low fields, the magnetic field 
dependence is rather similar to a purely diamagnetic case in which the density of states rises 
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smoothly and saturates at some critical value (at v = 1.) Again, this is softened, there is 
additional structure once the paramagnetism becomes comparable to the diamagnetism [i.e. 
a V.) For large enough fields the paramagnetism is always dominant. This means that 
N/ N will always approach unity from above. 

Unfortunately, most of this discussion is academic. As argued above, the cross-over field 
where a = v for if || c is about lOOOT. The value where a = 1 is about 350T while the value 
where i/ = 1 is about 70T. These are the scales where these affects would be visible, and 
they are probably beyond any experimental reach. For H \\ ah ii is reasonable to expect a 
similar structure to that shown in Fig. The field where a = z/ is about lOOT. The value 
where a = 1 is still 350T and the value where z/ = 1 is about lOOOT (again, depending on 
the value chosen for a). Again, all of these fields are beyond experimental reach and so the 
results of this appendix could probably not be observed in experiment. Nevertheless it is of 
theoretical interest to understand the analytic structure of the density of states. 
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